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DESCRIPTION ^^^^ 

Method of Optically Analyzing Inhomogeneous Medium 
Technical Field 

The present invention relates to a method of optically 
5 analyzing an inhomogeneous medium. 

Background Art 

While natural products, organisms, etc. typify 
inhomogeneous scattering media, techniques for analyzing such a 
scattering medium which scatters light while having pi nonuniform 

10 structure have not been known well. Of course, inhomogeneous media 

can approximately be analyzed as a uniform subject, though. 

According to Japanese Patent Application Laid-Open No. HEX 
2-234048, an apparatus has been developed, which irradiates an 
organism with high-frequency pulse light (on the order of 70 MHz 

15 to 200 MHz) having several discrete wavelengths (e.g., 760 nm, 

etc.), receives the light at a certain distance, determines 
absorption and scattering coefficients at a measurement wavelength 
by using a diffusion equation with respective amounts of the 
attenuation, phase deviation, and modulation of the light, and 

20 estimates theamount of hemoglobin from the absorption coefficient . 

Disclosure of the invention 

However, in an actual scattering medium, optical 
characteristics and refractive indexes may vary among its tissues, 
thereby causing mutual scattering among the tissues. Namely, 

25 biological tissues are not uniform, and cannot be obtainedby simply 

mixing scatterers and absorbers uniformly. In an organic tissue. 
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its characteristic substances (bloodvessels, muscles, fats, etc.) 
are distributed in a complicated fashion, so that optical 
characteristics {e*g., transmittance distribution) of the 
organism differ from those of its ideal uniform system. However, 
5 the conventional method calculates the optical characteristics 

while assuming the organism to be a uniform system, and thus is 
hard to measure concentrations of substances in the organism 
accurately. 

More specifically, the conventional theory cannot describe 
10 behaviors of light in a medium when optical characteristics vary 

among individual tissues so that mutual scattering occurs among 
the tissues. When the refractive index of the medium varies, its 
light velocity varies as well, whereby an analyzing method assuming 
that light propagates at a fixed speed does not hold. Analyses 
15 such asmicroscopicMonte Carlo method are not commonly applicable, 

since they cannot be conducted unless a nonuniform structure 
distribution is assumed with precision. Therefore, such a 
conventional method cannot analyze physical properties 
accurately - 

20 Japanese Patent Application Laid-open No. HEI 4-297854 

discloses a physical quantity analyzing technique using an 
exponential function, which corrects light-diffusing effects in 
media . This method uses the exponential function as an absorbance 
correcting function, whereby effective optical path length 

25 increases because of light diffusion even in a uniform medium. 

Optical intensity attenuates when the effective optical path length 
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increases , 

Therefore, the physical quantity analysis of Japanese Patent 
Application Laid-OpenNo. HEX 4-297854 empirically introduces an 
exponential function having a constant with an unclear physical 
5 content in order to modify the apparent optical intensity when 

correcting influences of changes in absorbance. 

For overcoming such a problem, it is an object of the present 
invention to provide a chemical analysis method which can 
accurately analyze inhomogeneous media* 

10 The inventor of the present application has studied, as a 

method of chemically analyzing an inhomogeneous medium, an analysis 
method which expresses behaviors of light in two tissues (x, y) 
with respect to a light advancing direction when setting up an 
initial differential equation concerning an inhomogeneous model. 

15 Such an analysis method is theorized according to an increase in 

scattering fed from another channel, while being based on the 
attenuation of light caused by absorbent substances and scattering 
substances . 

This differential equation and a conventional differential 
20 equation (e.g., Kubelka-Munk) basically differ from each other 

in their setup, so that the absorption and scattering coefficients 
in a general expression obtained by the foregoing usually differ 
from those in a theory used in general . It is therefore difficult 
for such a technique to evaluate results obtained by directly using 
25 a general expression in comparison with those obtained by the 

conventional method. 



3 



2004f 6^258 mm 80EI PATET i LAW FIRM NO. 8175 P. 17 

b'P0;^-U261-00 



The method of the present invention conforms to a 
Kubelka-Munk type analysis method, and is based on an equation 
with respect to a light advancing direction, two simultaneous 
equations with respect to an opposite direction, and a general 
5 expression solving a differential equation describing the giving 

and taking of photons between two phases (x, y) . 

Therefore, when the part of mutual giving and taking of 
photons between the phases concerningnonunif ormity is eliminated, 
a conventional Amy expression (Kubelka-Munk type general 
10 expression) is obtained, whereby the effectiveness of this 

nonuniformity theory can be proved more easily than the 
conventional method, and the most common solving method is yielded . 

Characteristic features of the present invention lie in that 
a general expression explaining nonuniformity is determined, and 
15 that the general expression can determine an absolute value in 

a target tissue even when a medium itself has a complicated 
nonuniform area in which tissues are intertwined • 
Best Modes for Carirying Out the Invention 

In the following, the method of chemically analyzing a 
20 inhomogeneous medium in accordance with an embodiment will be 

explained. Constituents identical to each other will be referred 
to with numerals or letters identical to each other without 
repeating their overlapping explanations. 

In this method, a solution of simultaneous equations 
25 describing a state receivingoptical scatteringbetweenaplurality 

of tissues is given in a function form in a method of chemically 
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analyzing an inhomogeneous medium, which makes light incident on 
the inhomogeneous medium, detects the intensity of light emitted 
therefrom, inputs thus detected optical intensity into a 
predetermined function, and determines the physical quantity of 
5 the inhomogeneous medium. An undepicted arithmetic unit 

calculates a specific physical quantity by inputting the detected 
optical intensity and a known coefficient in a database into the 
function. 

This function is defined by absorption coefficients (a) and 
10 scattering coefficients (s) of respective media of the tissues. 

Though this function also includes a tissue cross-sectional area 
ratio (K) and a mutual scattering coefficient (F) which are 
characteristics of nonunif ormity, they are not required to have 
explicitly been known, which is a characteristic feature of this 
IS method* 

In a specific example in which a medium is composed of two 
kinds of tissues (X, Y) , physical quantities concerning them are 
distinguishedby suffixes (x, y) . Namely, suf fixes x and y indicate 
that their absorption coefficients a and scattering coefficients 

20 s relate to the respective tissues X and Y. The numbers of photons 

traveling the tissues X and Y are represented by x and y. The 
light entering direction is defined as Z direction. Since the 
photons are also converted into an opposite flow by scattering 
during their progress, those traveling forward and backward are 

25 provided with suffixes f and b, respectively. The numbers of 

photons traveling forward are Xf and yf, whereas the numbers of 
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photons traveling backward are Xb and yb. 

Basic equations are given by the following: 
Behavior of light concerning the first tissue X; 



Behavior of light concerning the second tissue Y: 

The coefficient determinant for expressing these 
10 differential equations is composed of four rows by four columns • 

If there are three kinds of tissues, expressions are made in 
conformity thereto, so as to form a coefficient determinant of 
six rows by six columns by using three kinds of optical coefficients 
and three mutual light quantity redistribution coefficients. 
15 These expanding methods have been known. 

The above expressions are employed when the average radius 
of curvature r of a tissue interface is greater than the photon 
mean free path i, i.e., when r > i - When r < J, the above-mentioned 
FSx or FSy is simply replaced by F. 
20 When the above-mentioned differential equations are solved, 

transmittance (T) and reflectance (R) are given by the following 
expressions: 
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exp(- A^/z) + — exp(- A^h) 



r.4 



7£-[(l H- - (1 - a exp(- 2A,A)]^ ^ a,)^ - & exp(- 2A,a)] 

- )& - exp(- 2A,/»)] + (1 - gjl - exp(2 - A/)j 

^[(UaJ^-(l-aJexp(-2A^)]+^[(l+a,)^-(l-a,^xp(-2A,/,)] 



where Ax/ Ay, px, py^ Ox/ and Oy are as follows: 

Since the depth (h) is longer than the photon mean free path 
length, the above-mentioned expressions are simplified as follows : 

r-4 i — ^ — ^ 



Jt 



Its two solutions are defined as Ax and Ay. 
10 When 

the following expressions are effective by an approximation . 

Therefore, the present invention also encompasses the case 
where the following expressions are used. 



15 A. = 
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or « 

5 As mentioned above, F^SxSy is simply replaced by when r 

< i here as well. K is the cross-sectional area ratio of X and 
Y tissues, whereas F is the mutual scattering light exchange 
coefficient between tissues (light quantity redistribution 
coefficient by mutual scattering between X and Y tissues) . None 

10 of them is required to be given explicitly as a value as will be 

explained later. 

If all of a, s, K, and F for the tissues are known, T and 
R can be determined- In a measurement analysis, on the other hand, 
necessary values among a, s, K, and F can be determined if values 

15 of T and R are measured. Namely, since reflected light intensity 

IiR) = lo ^ R or transmitted light intensity I(t) = lo T can be 
measured, each of them becomes an expression with a, s, K, and 
F alone if incident light intensity lo at a known wavelength is 
known. If K and F are known as coefficients, for example, two 

20 simultaneous equations including unknowns a and s can be obtained. 

Unknown coefficients, i.e., target physical quantities, can be 
calculated if the known coefficients are stored in a database and 
measured values are put into the above-mentioned functions in the 
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arithmetic unit together with the known coefficients. 

In this case, the unknowns can be determined if equations 
whose number is not smaller than that of the unknowns are set up. 
While it has conventionally been impossible to calculate 
5 concentrations of ingredients from measured results in the case 

of inhomogeneous media, this method makes it possible to calculate 
accurate ingredient concentrations for the first time, 

F will now be explained. In general, in order for photons 
scattered within tissues on both sides of an interface therebetween 
10 to enter their neighboring tissues beyond the interface, their 

scattering points must exist within the photon mean free path from 
the interface. The volume concerned exists within the mean free 
path from the interface. 

As the scattering coefficient is higher, the quantity of 
15 light mingling into the neighboring tissue upon scattering 

increases . On the other hand, the mean free path becomes shorter, 
whereby the total volume of related generation points decreases. 
The increase and decrease cancel each other out, so that quantities 
of light mingling from neighboring tissues are equal to each other 
20 and are constant regardless of the scattering coefficient. 

This can be proved by mathematical calculations as well. 
Though the above relates to a case where the interface between 
tissues is a plane, a case where the interface between tissues 
is a curved surface will be considered separately therefrom. Here, 
25 it is assumed that one of the tissues has a curved surface such 

as that of a capillary in an organism, while its average radius 
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of curvature is shorter than the photon mean free path on the contrary 
to the above. 

The light scattered in a blood vessel always exits to an 
external tissue^ so that its related generation voliime is the blood 
5 vessel volume itself. On the other hand, generation points of 

photons which are scattered in tissues on the outside of the blood 
vessel and mingle into the blood vessel only exist on the outer 
surface of a cylinder surrounding the blood vessel while having 
a radius therefrom equal to the mean free path. 

10 The product of the related volume and the solid angle by 

which scattered light diverging therefrom just reaches the inside 
of the blood vessel becomes an effective volume for exchanging 
photons. When calculated, this volume can be proved to be 
sxibstantially equal to the original blood vessel volume. 

15 As explained In the foregoing, the above-mentioned 

differential equations are based on the inventor's finding that 
the mutual mingling coefficient of light between inhomogeneous 
media is simplified and can be considered constant. However, the 
equations explained in the foregoing relate to a case where the 

20 radius of curvature of the interface is shorter than the mean free 

path. When the radius of curvature of the interface is longer, 
by contrast, an expression in which F^SxSy in the above-mentioned 
equation is simply replaced by another constant is used. It 
will be sufficient if these equations are selectively used 

25 depending on the object. 

The tissue cross-sectional area ratio K is originally a 
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geometric constant and can be considered an undetermined constant . 

However, there may be a case where the ratio is not constant 
in the light penetrating direction z in a medixim to be measured. 

In such a case, the tissue cross-sectional area ratio K is 
treated as an average constant concerning the depth of the medium . 
When K changes moderately with respect to the depth, this change 
can be neglected. 

The number of equations can be increased by changing 
measurement wavelengths. Of course, the wavelength dependence 
of each of coefficients in each case is required to be known. 

The foregoing technique can determine the optical absorption 
coefficient a, vertical scattering coefficient s, etc. 

In general, physical quantities to be determined by a 
measurement analysis are ingredient concentrations and the like. 
Here, if two ingredients {p, q) exist in both of the tissues X 
and Y, the absorption coefficient a and scattering coefficient 
p are given by the following expressions (*) : 
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Since a/p, a/q, s/p, s/q, and the like are absorption or 
scattering coefficients per unit concentration of their 
corresponding target ingredients (in the mediiim) and are known 
in general, p and q can be determined by an algebraic calculation 
5 from a and s determined as above* 

Since constants can be erased later, the measurement 
background noise (constant) may be incorporated in addition to 
the above-mentioned p and q proportional terms when necessary. 
These equations are based on the proportional addition rule 
10 of the most common concentration media. 

In practice, scattering coefficients are likely to change 
in high concentration media because of multiple scattering or 
association of scattered particles and the like, so that the 
proportionality may be lost. 
15 When multiple scattering occurs, the total scattering tends 

to be saturated. At this time, a proportionality such as that 
of the following equation, which is weaker than the linear 
proportionality, is employed, for example- Which proportionality 
to employ is determined beforehand by a study depending on the 
20 object. 

or 5=-^ 

Meanwhile, the above-mentioned expressions include unknown 
coefficients (K, F) . These coefficients are known to be constants 
and thus can be erased easily. It will be sufficient if measurement 
25 is per formed with varied wavelengths in order to set up an additional 
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number of equations corresponding to the number of un3cnown 
coefficients. 

In general/ however, there is a high risk of increasing errors 
if the minimal niomber of equations are set up while relying on 
5 the wavelength dependence alone. This is because errors often 

exist in various coefficients themselves and accumulate in final 
results. It is therefore important to carry out measurement in 
many wavelengths other than those minimally employed, in order 
to alleviate the errors. In this case, arithmetic operations are 

10 performed according to the method of least squares. 

When determination equations are nonlinear as in the present 
invention, least-squares calculations of nonlinear equations are 
carried out according to Newton's approximation. Here, when 
arithmetic operations are complicated as in the present invention, 

15 calculations in which the least-squares method is applied to 

individual stages for respective equations mentioned above seem 
to be easier for determining final p and q. 

When calculations are employed in a divided fashion as such, 
the convergence of errors, which is a characteristic feature of 

20 the least-squares method, is divided into the respective stages, 

whereby there is no guarantee that an optimal value is finally 
obtained. 

In this regard, since all the stages are expressed by 
analytical expression? though nonlinear, the present invention 
25 can analytically express the coefficient matrix of the determinant 

in each stage and can attain the total solution by simply applying 
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10 



15 



the least-squares method to the final single determinant, while 
being able to derive coefficient determinants from the individual 
stages and combine them together, which are excellent features . 
In the following, an example of solving method will be set forth • 
Since Newton's approximation is well known, only its 
fundamental means will be explained in brief. When (e.g,, 16) 
values measured with respective wavelengths different from each 
other of deviations from an estimated resolution value taken as 
a center are referred to as ml, ml6, they are represented 

as in the following equation since T and R are described by p and 
q, [M] is a Jacobian determinant according to differential classes 
of the individual coefficients and is a known coefficient. 



^16 
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[M] is diagonalized when necessary. 

Siobsequently, matrix [A] is determined by the 
above-mentioned expression (*) , whereas a and s are described by 
p and q, whereby the following equation is obtained similarly: 
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Hence, the solution of the following equation can be 
obtained: 



Anil 
Am 
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In these numeric calculations, accelerated convergence 
5 solving techniques can be employed by various known methods. In 

the case of the above-mentioned example, calculations are possible 
withanynumber of wavelengths (e.g., 4096 wavelengths) not smaller 
than 6, regardless of whether they have a continuous spectrum or 
discrete wavelength characteristics. If a sufficient number of 
10 measurement operations are carriedout with respective wavelengths 

different from each other, final calculation results can be 
obtained with the least error, which has conventionally been 
impossible in such a complicated system^ 

In particular, there has conventionally been no theory which 
15 enables nxomeric calculations by such an analytic technique. It 

is a new discovery that the interrelationship between nonuniform 
minute tissues can be expressed and processed by simple 
coefficients of K and F. 

If it is clear that the medium has totally different layer 
20 structures in its depth, e.g., where a skin layer, a coating 

protecting layer, and the like exist, they are preferably divided 
for respective calculations. The calculations are simple when 
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the partial thickness of each layer is known, but can also attain 
solutions when the thickness is unknown. Namely, the coefficient 
matrices explained above can be combined together, so as to be 
treated as a cascade of two kinds of transfer functions. 
5 It is indicated that calculations should include some 

undetermined coefficients in the case explained in the foregoing. 
It is shown that complicated tissue structures in inhomogeneous 
media and the like can be expressed as a form including unknown 
coefficients by setting up a number of simultaneous equations and 

10 solving them by employing the least-squares method. A 

characteristic feature of this theory also lies in that solutions 
can thus be determined while utilizing the fact that unknown 
coefficients can be erased as implicit functions. 

Therefore, metrical solutions can be determined even when 

15 a tissue has a nonuniform structure and a complicated, unknown, 

intricate form or an unknown cross-sectional area. 

This immediately guarantees quantitative measurement 
results of p, q, and the like. For example, though conventional 
biological optical analyses can easily detect the absorption of 

20 oxygenated hemoglobin in the blood, they cannot determine an 

absolute ingredient amount quantitatively but measure only 
qualitative behaviors if the effective cross-sectional area ratio 
of the blood vessel is unknown. By contrast, the method of the 
present invention makes it possible to erase the effective 

25 cross-sectional area and the like and determine the absolute value 

of a specific ingredient in a specific tissue for the first time. 
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Thus not only actual objects/ organisms ^ and the like are 
nonuniform, but also their tissues are complicated and have 
conventionally been impossible to display and analyze 
mathematically. The technique of the present invention has found 
5 that tissue structure parameters can be approximated as simple 

implicit functions, and erases them in the process of a nxzmeric 
analysis/ thus exhibiting a characteristic feature of being able 
to reach solutions without arguing nonuniform contents in detail. 
When a substance is a medium having a nonuniform structure 

10 or a multiphase or multiphase separation structure as in a mixture, 

the above-mentionedmethodcanbe applied thereto, so as to separate 
individual components in each tissue from each other and calculate 
them independently from each other . When the inhomogeneous medium 
is composed o'f a number of sectional tissues, simultaneous 

15 equations including the number of intermediate elements 

corresponding to the kinds of tissues can be contracted, and thus 
contracted equations can be applied to the above-mentioned 
analyzing technique in the case where the number of tissues X, 
Y is 2. 

20 The analyzing method explained here is quite useful in 

practice, since target solutions canbe attainedwithout explicitly 
handling differences in detailed optical characteristics among 
a plurality of tissues composed of tissues X, Y, etc. In a 
measurement example in an organism as a specific example, blood 

25 vessels and other tissues exist as an inhomogeneous multiphase 

mixture state in an actual organism. According to the example 



17 



20Hf 6^258 1 5^45i) SOEI PATET & LAW FIRM 



i;T02-u261- 



mentioned above, the X tissue is treated as the blood part in the 
blood vessels, the Y tissue is treated as* the other tissue part, 
andx and y are calculated as respective quantities of light included 
therein. 

5 Measurement of concentrations of oxygenated hemoglobin and 

reduced hemoglobin in the X tissue, measurement of concentrations 
of protein groups (Cytochrome, aa3, etc.) of electron transport 
systems in Y tissue cells, etc. become possible . 

Such characteristic features make the present invention 
10 effectively applicable to media which are mixed nonunif ormly with 

absorption and scattering coexisting, and a wide range of 
substances conventionally hard to measure, thus yielding great 
effects. Examples of the substances include industrial products, 
natural products such as food and plants, other mixtures, 
15 polycrystals, and turbid materials. 

In brief , the present invention is effective when substances 
having respective modes different from each other are mixed 
nonun i f o rml y . 

As for foods, while cereals, fruits, oils and fats, etc. 

20 are vigorously studied and developed in terms of quantitative 

measurement of inner substances and the like by utilizing infrared 
light, analyses are in progress while their samples are considered 
uniform. However, foods include those which should clearly be 
treated as nonuniform substances, such as a product in which cut 

25 fruits are contained in a yogurt. The method of the present 

invention wor.ks effectively for them. 
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In the field of color matching of fiber dying (color matching 
calculations) , the method of the present invention is not effective 
in the field of color matching calculations for a scheme in which 
fabrics are dyed after being prepared, but in the field of so-called 
5 fiber blend color matching calculations in which strings or wool 

yarns are dyed beforehand and then these strings having respective 
colors slightly different from each other are mixed or strings 
having respective colors different from each other are mixed so 
as to yield a desirable color. 

10 In the field of water, plastics, and other (metal) materials 

(while assuming visible light to be a common electromagnetic wave 
instead of medium) , a solid, a liquid, or their transient state 
(solid solution in the case of a metal) exists in their melted 
stage. The present invention can become means for elucidating 

15 details of their mixed state. 

Also, the mathematical theory solely extracted from this 
specification is seen to be effective in the study and elucidation 
of nonuniform and inhomogeneous systems in future. This is because 
the science of nonuniform systems has just begun at present, so 

20 that substantially no theoretical expression for explaining a 

nonuniform system exists except for that of the present invention. 

For example, when a developing process in which an ovum grows 
to a human by way of an embryo is optically studied, nonuniformity 
increases . 

25 The present invention canbecome one of mathematical theories 

which gives a hint to a grand problem of why an organism makes 
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a nonuniform system and complex system so skillfully and perfectly, 
and an experimental system therefor in future. 

However, it is necessary that (i) a numerical database of 
absorption and scattering coefficients be known and that (ii) the 
5 convenience of numerical processing be studied since the method 
is based on a generic solution, whereby the amount of calculations 
can increase. These may vary depending on objects, whereby it 
is necessary to find appropriate simplifying means empirically 
for specific cases. 
10 While it is stated that the parameters of the tissues can 

be approximated as simple implicit functions and erased in the 
process of numerical analysis, it is one of points in practice 
how calculations are carried out without letting out the tissue 
parameters. The above-mentioned arithmetic operations are 
15 performed by a computer as an arithmetic unit. 

The expression on which the above-mentioned general 
expressions are based is an expression of Amy (L.Amy, Rev. d'optique 
16 81 (1938) ) . The general expressions of the present invention 
are based on Amy' s expression (uniform composition expression for 
20 a single composition) and developed into a nonuniform state, i.e», 

giving and taking of photons between different channels are added 
thereto. Thus designed expressions can be used as general 
expressions for quantitatively analyzing substances in mixtures, 
polycrystals, turbid materials, and other nonuniform states. 
25 As explained in the foregoing, the method of optically 

analyzing an inhomogeneous medium in accordance with the 
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above-mentioned embodiment determines a physical quantity of a 
light-absorbing and light-scattering inhomogeneous medium 
including a plurality of tissues X, Y having respective optical 
characteristics different from each other from the intensity of 
5 light transmitted therethrough or reflected thereby. Aprocedure 

executed by a computer comprises the following steps: 

(i) Light having a known intensity lo at a known wavelength is 
made incident on the medium from a light source. 

(ii) The optical intensity I{r, upon reflection or optical 
10 intensity I{t) upon transmission emitted from the medium is measured 

by a photodetector . 

(iii) The above-mentioned steps (i) and (ii) are executed at a 
plurality of wavelengths, and the optical intensity I(R) or I(t) 
as a measured result is recorded into a storage device of the 

15 computer. 

(iv) Using a known physical quantity (PV) and the optical 
intensity I(r) or I(t) recorded in the storage device, a target 
physical quantity is calculated by an arithmetic operation of the 
computer . The tissues are composedof two phases of X and Y, whereas 

20 the target physical quantity is set to a physical quantity other 

than the known physical quantity (PV) in the database. 

Namely, the target physical quantity is determined by 
inputting those known in the following coefficients as known 
phys ical quantities stored in the database and the measured optical 

25 intensity I (r^ or I (tj into those related to the arithmetic operation 

in the following expressions: 

21 

I 
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(1) detected value I(r) or I(T) 

(2) physical quantities in the database 
T: transmittance of the medium 

R: reflectance of the medium 
5 h: thickness of the medium 

ax: absorption coefficient of tissue X 
Sx! scattering coefficient of tissue X 
ay: absorption coefficient of tissue Y 
Syi scattering coefficient of tissue Y 
10 pxt parameter defined for tissue X in conformity to a 

uniform system 

py: parameter defined for tissue Y in conformity to a 
uniform system 

Ox: pa=rameter defined for tissue X in conformity to a 
15 uniform system 

Oy: parameter defined for tissue Y in conformity to a 
uniform system 

K: cross-sectional area ratio of tissues X and Y 
F: light quantity redistribution coefficient of mutual 
20 scattering between tissues X and Y 

Ax/ Ay: two solutions of A equation. 

Here, each of (axf Sx) and (ay, Sy) has conventionally been 
written as (pa/ Ps' ) in academic literatures. 

(3) Expressions 

25 reflected light intensity I(r) = To x R 

transmitted light intensity 1(t) - lo ^ T 
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R 



X(l-crJ+(l-orJ 



T 



4 



* y 



X(l + aJ+(l+aJ 
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The mutual redistribution coefficient F is determined by 
fine structures of tissues and mostly is a constant. Namely, F 
in the above-mentioned expression is used when the average radius 
of curvature of the tissue interface is smaller than the photon 
mean free path, whereas another constant F2 is simply used in place 
of F2SxSy in the expression when the average radius of curvature 
is greater than the photon mean free path. 

Using relationships in which the physical quantities ax, 
Sxf ay, and Sy are proportional to the concentration of a 
predetermined target ingredient in the medium, the concentration 
of the predetermined ingredient can be calculated. 

When a nonproportional relationship due to multiple 
scattering or the like exists between the concentration of each 
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ingredient and its absorption or scattering light quantity instead 
of a usual proportional relationship, solutions can be determined 
by employing their relational equation which is not linearly 
proportional . 

5 The optical analysis method of the present invention can 

accurately analyze inhomogeneous media. 
Industrial Applicability 

The present invention can be utilized in a method of optically 
analyzing inhomogeneous media. 
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